REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  the  Department  of  Defense,  Executive  Services  and  Communications  Directorate  (0704-0188).  Respondents  should  be  aware 
that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB 
control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 


2.  REPORT  TYPE 

Conference  Proceeding 


1.  REPORT  DATE  (DD-MM-YYYY) 

29-04-2015 


4.  TITLE  AND  SUBTITLE 

EO  Signal  Propagation  in  a  Controlled  Underwater  Environment 


3.  DATES  COVERED  (From  -  To) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


6.  AUTHOR(S) 

Weilin  Hou,  Silvia  Matt 


5c.  PROGRAM  ELEMENT  NUMBER 


0602782N 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

Oceanography  Division 

Stennis  Space  Center,  MS  39529-5004 


73-6604-04-5 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NRL/PP/7 330-1 4-2 1 03 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 

One  Liberty  Center 

875  North  Randolph  Street,  Suite  1425 

Arlington,  VA  22203-1995 


10.  SPONSOR/MONITOR  S  ACRONYM(S) 


11.  SPONSOR/MONITOR  S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  is  unlimited. 


14.  ABSTRACT 

Underwater  electro-optical,  or  EO,  transmission  is  a  function  of  medium  properties  and  constituents  within.  While  the  majority  of 
the  research  focus  has  been  on  the  constituents,  especially  the  particulate  forms,  recent  research  indicates  that  under  certain 
conditions,  the  apparent  signal  degradation  could  also  be  caused  by  variations  of  the  index  of  refraction  associated  with  temperature 
and  salinity  micro  structure  in  oceans  and  lakes.  These  would  inherently  affect  optical  signal  transmission  underwater,  which  is 
important  to  both  civilian  and  military  applications  involving  search  and  rescue,  intelligence,  surveillance  and  reconnaissance 
applications,  as  well  as  optical  communications. 


15.  SUBJECT  TERMS 

Controlled  turbulence;  simulation;  Rayleigh-Benard  tank;  phase  screen;  underwater;  turbulence;  optical  scattering 


16.  SECURITY  CLASSIFICATION  OF: 
a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE 

Unclassified  Unclassified  Unclassified 


17.  LIMITATION  OF  18.  NUMBER  19a.  NAME  OF  RESPONSIBLE  PERSON 
ABSTRACT  OF  Weilin  Hou 

PAGES  - 

UU  .  19b.  TELEPHONE  NUMBER  (Include  area  code) 

6  (228)  688-5257 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


PUBUCATION-OR-RBESENTATION  RELEASE  REQUEST 


-  I'Z&drQLf'H 


NRLINST  5600.2 


Ref:  (a)  NRL  Instruction  6000.2 
(b)  NRL  Instruction  5510.40D 

End:  (1)  Two  oopte*  of  subject  paper 
(or  abstract) 


t  )  Abstract  only,  publfihed 
(  )Book 

(  )  Conference  Proceedings 

(refereed) 

!)  Invited  speai^r 

)  Journal  article  (refereed) 

)  Oral  Presentation,  published 
(  )  Other,  explain 


(  )  Abstract  only,  not  published 

(  )  Book  chapter 

( X  )  Conference  Proceedings 
(not  refereed) 

(  )  Multimedia  report 

(  )  Journal  article  (not  refereed) 

(  }  Oral  Presentation,  rot  published 


Route  Sheet  No  . 

Job  Order  No.  73*  6004-  04-  S  ^ 
aassIHcalion^* — -X  u 

Sponsor /XiQNR  BASe)(,?  , 

appmvaTobfekied^^  yes  X  no 


Title  of  Paper  or  Presentation 

£0  Signal  Propagation  in  a  Controlled  Underwater  Environment 


Authors)  Namejd)  (FfcsLAf/.Laad,  Cotef,  AftiHation  if  not  htRL 
WeUtn  Hon  7333  Silvia  Matt  NR**/ 


t  is  intended  to  offer  this  paper  to  the  MTS/IEEE  Queans  2014 . . . . . . . . _ 

(Name  cf  Conference) 

)7- APR  -  10*  APR- 14,  Taipei,  Taiwan,  Unclassified 


(DaHi,  Piece  and  Classification  of  Conference) 

and/or  for  publication  In  MTS4ESE  Oceans  2014,  Unclassified 

(Name  and  Classlfcatton  of  Publication)  (Name  of  Publisher) 

After  presentation  or  publication,  pertinent  publicationfpreeentatlon  data  will  be  entered  in  the  publications  data  base,  In  accordance 
rvlth  reference  (a). 

it  is  the  opinion  of  the  author  that  the  subject  paper  (is _ )  (is  not _ *)  classified,  in  accordance  with  reference  (b). 

This  paper  does  not  violate  any  disclosure  of  trade  secrete  or  suggestions  of  outside  individuals  or  concerns  which  "have  been 

communicated  to  the  Laboratory  in  confidence.  This  oaper  (does _ )  (does  not _ X)  contain  any  militarily  criticaj^technology. 

This  subject  paper  (has  _ _ )  (has  never  X  )  been  incorporated  in  an  official  NRL  Report 

„ _ WelUn  Hou.  7  ^ 


Name  and  Code  (Principal  Author)  (SJg 


mraari 


Meed  toy  ^  / d 


PubHdy  ecceasitote  soared  used  brlhla  puWtealion 


3  ■'J&3J 'jKsV 


C  cJv>  ") 


ADOR/Oirector  NCST 

6.  R.  French!,  7000 


Public  Affairs  (Unctexifod/ 
Unlhnitad  Only),  Code  7030  4 


Division.  Code 


Author,  Code 


IQ-NRl  *511*  (Rev.  12-W)  fe) 


Paper  or  afc&art  wag  <efeas«d 
A  copy  is  filed  In  this  office, 


nullify  the  Securtty^tev 


THIS  FORM  CANCELS  AND  SUPERSEDES  ALL  PREVIOUS  VERSIONS 


t>0/£0  3Wd 


9h06 L3LZQZ 


Lb'Ll  t>T0Z/0I/E0 


EO  Signal  Propagation  in  a  Simulated  Underwater 

Turbulence  Environment 


Weilin  Hou 

Naval  Research  Laboratory 
Hydro  Optics,  Sensors  and  System  Section 
Stennis  Space  Center,  MS,  USA 

Abstract —  Underwater  electro-optical,  or  EO,  transmission  is  a 
function  of  medium  properties  and  constituents  within.  While  the 
majority  of  the  research  focus  has  been  on  the  constituents, 
especially  the  particulate  forms,  recent  research  indicates  that 
under  certain  conditions,  the  apparent  signal  degradation  could 
also  be  caused  by  variations  of  the  index  of  refraction  associated 
with  temperature  and  salinity  microstructure  in  oceans  and 
lakes.  These  would  inherently  affect  optical  signal  transmission 
underwater,  which  is  important  to  both  civilian  and  military 
applications  involving  search  and  rescue,  intelligence, 
surveillance  and  reconnaissance  applications,  as  well  as  optical 
communications. 

To  study  the  effect  of  optical  turbulence  and  to  mitigate  its 
impacts,  a  controlled  environment  allowing  various  intensities  of 
turbulent  mixing  is  a  critical  asset.  Numerical  experiments  as 
well  as  measurements  have  been  carried  out  in  such  a  simulated 
environment,  in  order  to  understand  mixing  setup  time, 
development  and  dissipation  rates.  The  domain  is  modeled  after  a 
large  Rayleigh-Benard  convective  tank  with  a  length,  width  and 
depth  dimension  of  5,  0.5  and  0.5m,  respectively.  The  convective 
mixing  is  realized  by  using  heating  and  cooling  plates  at  the 
bottom  and  top  of  the  tank  at  given  temperature  differences.  The 
computational  fluid  dynamics  model  is  implemented  with  large 
eddy  simulation  approximation.  Dissipation  rates  from  model 
and  measurements  are  compared  and  suggest  fully  developed 
turbulence  has  been  achieved  by  this  setup.  Optical  signal 
transmission  under  these  conditions  are  also  examined,  through 
image  degradation  using  image  quality  metric,  and  phase  screen 
models  from  corresponding  power  spectrum.  The  integrated 
temperature  variation  along  the  transmission  path  is  compared 
to  generated  phase  screens,  along  with  discussions  on  reducing 
uncertainties  in  estimation  of  key  parameters. 

Keywords-component;  controlled  turbulence;  simulation ; 
Rayleigh-Benard  tank;  phase  screen;  underwater;  turbulence; 
optical  scattering 

I.  Introduction 

Recent  research  on  underwater  vision  and  optical,  as  well  as 
acoustical  signal  propagation  suggests  better  understanding  is 
needed  to  quantify  or  mitigate  the  influence  of  underwater 
optical  turbulence  [1-5].  By  its  definition,  turbulence  is  chaotic 
in  nature,  where  statistical  approach  has  been  the  only 
effective  tool  in  quantifying  related  processes,  along  with 
idealized  conditions  and  assumptions.  Observations  have  been 
made  both  in  lab  and  field,  along  with  ancillary  measurements, 
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in  order  to  have  a  reasonable  estimate  of  the  processes 
involved.  The  majority  of  the  turbulence  research  has  been 
carried  out  in  the  atmosphere,  for  studies  related  to  weather 
events  in  general.  Specific  to  optical  signal  transmission,  the 
efforts  are  largely  stemmed  from  the  need  to  acquire  better 
astronomical  images  of  distant  stars,  which  prompted  the 
seminal  paper  by  Fried[6],  that  the  seeing  is  dependent  of 
optical  turbulence  intensity,  where  it  can  be  represented  by  the 
Fried  seeing  parameter  r0.  If  the  wavelength  of  light  passing 
through  the  system  is  X,  then  the  optical  transfer  function, 
OTF,  or  modulation  transfer  function  when  only  amplitude  is 
involved,  is 

5 

HQV)  =  exp  [-3.44  (^f]  (1) 

here  W  is  the  spatial  frequency  in  cycles  per  radian,  The 
exponential  term  in  the  OTF  is  essentially  the  structure 
function  of  the  phase  variation,  due  to  the  index  of  refraction 
change  caused  by  the  temperature  variation  in  the  atmosphere. 
This  relates  back  to  the  basic  statistical  imaging  principle, 
where  the  interference  between  two  field  points  forms  the 
image  that  describes  the  spatial  coherency.  This  is  known  as 
the  Van  Cittert-Zernike  theorem  [7,  8],  and  essentially  the  two 
dimensional  view  of  the  famous  Young’s  interference 
experiment.  The  spatial  perturbations  during  the  wave 
propagation  alter  the  imaging  outcome  at  the  pupil  plane,  by 
altering  the  interference  fringe  patterns.  Due  to  the  nature  of 
turbulent  flow  in  the  atmosphere,  where  the  motion  follows 
lateral  wind  directions  and  density  structures,  a  layered 
approach  can  be  used  for  most  astronomy,  as  well  as 
reconnaissance  imaging  needs.  However,  this  is  not  the  case 
for  the  most  oceanic  turbulence  study,  unless  vertical  optical 
transmissions  through  strong  horizontal  layers  are  involved. 

A  few  assumptions  have  to  be  made  to  properly  describe 
the  seemingly  random,  chaotic  process  loosely  termed 
turbulence.  Even  from  the  observations  from  early  days,  like 
those  of  Da  Vinci  and  his  water  swirl  drawing,  we  understand 
that  turbulence  relates  to  the  way  energy  dissipates  across 
various  spatial  scales.  Kolmogorov’s  classical  theory  of 
turbulence  is  often  used,  where  fully  developed  isotropic 
homogeneous  turbulence  conditions  are  assumed,  and  the 
energy  spectrum  follows  the  -5/3  law  in  the  inertial  sub-range, 
based  on  dimensional  analysis.  This  is  the  condition  for  which 


U.S.  Government  work  not  protected  by  U.S.  copyright 


the  above  equation  was  derived,  and  we  can  clearly  see  the 
trademark  signature  in  the  equation.  While  sounding 
somewhat  contradictory,  a  wide-sense  stationary  condition  has 
to  be  reached,  in  a  statistical  sense  over  temporal  (stable)  and 
spatial  scales  (isotropic),  in  order  to  describe  the  kinetic 
energy  dissipation,  as  well  as  other  scalar  quantities. 

Until  recently,  the  standing  view  about  the  source  of 
underwater  electro-optical  (EO)  imaging  degradation  has  been 
the  scattering  by  the  medium  itself  and  the  constituents  within, 
namely  particles  of  various  origins  and  sizes.  Recent  research 
indicates  that  under  certain  conditions,  the  apparent 
degradation  could  also  be  caused  by  variations  of  the  index  of 
refraction  associated  with  temperature  and  salinity  micro¬ 
structures  in  oceans  and  lakes  [1,  5,  9].  These  would 
inherently  affect  optical  signal  transmissions  underwater.  Two 
of  the  key  parameters  commonly  used  to  describe  turbulence 
are  directly  related  to  the  optical  turbulence  intensity,  S„,  [1, 
5],  a  coefficient  directly  related  to  the  optical  transfer  function 
(OTF)  of  the  turbulence  impact  (see  Eq.  3  below),  shown  in 
the  form  below,  including  path  radiance  and  particle  scattering 
contributions: 

OTF(r,r\otal=OTF(iy,r)^ 


here  90  relates  to  the  mean  scattering  angle,  c  and  b  are  the 
beam  attenuation  and  scattering  coefficients  respectively,  r  is 
the  imaging  range,  and  D  relates  to  path  radiance  [3].  Sn 
contains  parameters  that  are  dependent  on  the  structure 
function,  which  can  be  further  expressed  in  terms  of  the 
turbulence  dissipation  rate  of  temperature,  salinity  and  kinetic 
energy,  assuming  the  Kolmogorov  power  spectrum  type: 

Sn~ze~m  (3) 

where  s,  %  represent  the  kinetic  energy  and  temperature 
variance  dissipation  rates,  respectively.  These  are  the  key 
parameters  to  be  determined  from  numerical  simulations,  as 
well  as  measurements. 

Only  a  simple  case  of  a  point  source  response  to  the 
imaging  system  needs  to  be  examined,  as  all  image  formation 
can  be  treated  as  a  convolution  between  a  point  source  and 
system  transfer  function  [1,  10].  For  continuous  wave  front 
distortions,  as  a  result  of  turbulence  degradation  satisfying  Eq. 
1,  it  has  been  shown  that  a  Fourier  transform  pair  links  the 
phase  fluctuation  and  corresponding  power  spectrum  [7,  8]. 
One  can  then  sample  the  power  spectrum  (in  this  case, 
temperature  variance  dissipation)  to  create  phase  screens, 
which  represents  layered  distortions  or  integrated  effects.  The 
covariance  approach  will  not  be  discussed  in  this  initial  effort. 


To  better  understand  turbulent  micro  structure,  especially 
those  associated  with  temperature  and  salinity  variations  in  the 
underwater  environment,  it  would  be  ideal  if  a  repeatable, 
somewhat  controllable  environment  can  be  implemented.  This 
is  the  motivation  of  our  Rayleigh-Benard  (RB)  convective 
tank  setup  at  the  U.S.  Naval  Research  Laboratory  at  the 
Stennis  Space  Center,  MS.  We  will  first  discuss  the  setup  and 
lab  experiment,  and  numerical  simulations,  as  our  attempt  to 
establish  a  controlled  turbulence  environment.  This  will  be 
followed  by  observations  on  image  quality  degradation 
through  the  controlled  environment,  and  related  phase 
variations  under  different  power  distribution  schemes 
associated  with  the  controlled  environment. 

II.  Simulation  of  Controlled  Environment 

A  5  by  0.5  by  0.5  m  (L,  H,  W)  plexiglass  tank  has  been 
fitted  with  stainless  steel  plates,  in  order  to  control  the  surface 
and  bottom  temperature  for  generation  and  maintenance  of 
convective  mixing  (Figure  1).  This  classical  Rayleigh-Benard 
tank  allows  generation  of  optical  turbulence  at  various 
intensities,  shown  by  an  example  in  Figure  2.  The  left  side  of 
the  panel  displays  a  turbulence-free  image  through  filtered 
water  over  5m  in  pathlength,  under  a  beam  attenuation  value 
less  than  0.1m"1  measured  by  an  ac-9  (WETLabs).  Without 
particle’s  presence,  when  the  optical  turbulence  intensity 
increases  due  to  an  increase  of  the  temperature  gradient 
between  the  plates,  one  can  see  the  severe  degradation  of 
imaging  resolution  by  examining  the  resolution  chart  shown 
by  the  frame  on  the  right  of  Figure  2,  which  is  in  line  with 
results  of  Eq.2.  Note  that  in  the  foreground  of  the  image  is  a 
partial  image  of  an  Acoustical  Doppler  Velocimeter  (ADV), 
made  by  Nortek,  which  can  also  been  seen  in  Figure  1  (left). 
These  are  used  to  measure  the  small-scale  velocities  in  the 
tank  during  the  experiment,  along  with  fast  conductivity  and 
temperature  probe  (CT)  from  PME,  and  three  smaller  profiling 
ADVs,  Vectrino  (Fig.l,  right),  also  from  Nortek.  For  initial 
examination,  seeding  necessary  for  ADV  measurements  are 
added  to  the  tank  after  optical  measurement  using  the  same 
temperature  gradient  setup,  to  ensure  only  turbulence  effects 
on  optical  signal  transmission  are  in  place.  This  does  not  affect 
the  circulation  pattern,  once  the  gradient  is  set.  This  is 
confirmed  by  numerical  simulations  to  be  discussed  below. 
These  measurements  are  necessary  to  help  model  setup,  as 
well  as  validation. 


Figure  1.  Rayleigh-Benard  convective  tank  setup  at  the  U.S.  Naval 
Research  Laboratory,  Stennis  Space  Center.  Different  temperature 


gradients  are  implemented  by  two  stainless  steel  plates  on  top  and 
bottom.  Left  panel:  along  tank  view  of  target  through  filtered  water 
with  ADV  embedded;  Right  panel:  close-up  vie  /  of  CT  and  ADV 
setup  (see  text  for  details). 

The  governing  equation  of  any  fluid  is  straight-forward 
Navier-Stokes  equation,  shown  below  fo  ■  incompressible 
fluid,  in  the  absence  of  external  forcing  [11]: 

^  +  vw  =  -  —  +  vF2V  (4) 

dt  p  w 

where  v  is  the  kinematic  viscosity.  In  most  cases  it  is  not 
possible  to  arrive  at  an  analytical  solution  to  recover  the  flow 
field.  Numerical  solution  using  high  speed  computers  has 
proven  to  be  an  effective  tool.  The  domain  )f  interest  can  be 
subdivided  into  grid  points,  at  which  a  solution  can  be 
obtained  by  discretizing  the  governing  equations  (including 
Eq.  4,  in  addition  to  continuity  and  energy)  i  l  space  and  time, 
for  a  given  set  of  boundary  and  initial  conditi  >ns. 


Figure  2.  Sample  images  of  impact  of  optical  t  irbulence  over  5m 
pathlength  in  water,  viewing  resolution  chart  and  ADV  in  the 
foreground.  Left:  turbulence  free;  Right:  strong  optical  turbulence. 

However,  sacrifices  have  to  be  made  d  le  to  the  cost  of 
computation,  especially  at  fine  spatial  and  te  lporal  resolution, 
which  is  necessary  for  micro  structure  studies.  Scale 
discrepancies  between  different  physical  processes  demand 
approximations  to  be  made,  as  a  compromise,  even  with 
powerful  supercomputers.  For  example,  the  ptical  wave  field 
phenomenon  requires  sub-picosecond  resol  ition  to  be  fully 
resolved,  which  is  not  possible  in  practical  underwater 
applications  at  the  moment.  Even  at  this  rate,  it  is  not 
sufficient  to  include  the  spatial  resolution  requirement 
comparable  to  the  wavelength  of  light  (whic  ti  is  sub-micron), 
in  order  to  study  processes  involving  wave  field  interactions 
for  typical  visible  wavelength.  Another  exa  pie  can  be  found 
in  our  computational  fluid  dynamics  (CFD)  approach,  where 
the  scales  of  turbulent  flow  cover  several  orders  of  spatial 
magnitude  that  can  be  prohibitively  expensi  re,  to  fully  solve 
using  direct  numerical  simulation  (DNS)  of  Eq.  4. 
Traditionally,  several  parameterizations  hav  been  developed 
to  model  the  sub-grid  scale  turbulence  processes,  including 
Reynolds-averaged  Navier-Stokes  (RANS),  and  large  eddy 
simulation  (LES).  LES  is  essentially  a  1  >w-pass  filtering 
approach,  where  larger  scales  are  resolved  while  the  finer 
scales  are  modeled.  LES  explicitly  resolves  the  larger  eddies 
in  the  flow,  and  models  processes  small  :r  than  the  grid 
spacing.  This  improves  the  physical  representation  of  the 
turbulent  field,  and  is  adopted  for  our  work.  An  open-source 


CFD  package,  OpenFOAM,  i ;  used,  with  traditional 
Smagorinsky  model  as  a  sub  grid  scale  model.  Two 
resolutions  have  been  investigated,  on  at  1cm  (x,y,z)  and  one 
at  2.5  mm  (for  horizontal  x,  y)  md  5mm  (for  vertical,  z). 
These  are  equivalent  of  1.25  and  20  million  grid  points, 
respectively.  The  temperature  distribution  of  simulation  under 
fully-developed  (time=  1050s),  strong  turbulence  conditions  is 
shown  in  Figure  3.  The  index  of  re  faction  (IOR)  in  the  tank  is 
calculated  based  on  relationship  to  :he  temperature  fields  [12]. 
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Figure  3.  LES  simulation  temperature  (in  Kelvin),  vertical  velocity 
and  IOR  in  the  Rayleigh-Benard  convective  tank  under  strong 
turbulence,  when  fully  developed  (see  text  for  details). 

One  may  notice  that  several  1  rger  scale  convective  cells 
(~1.5m)  are  present  in  the  above  iimulation  output,  which  is 
confirmed  by  visual  observation;,  as  well  as  temperature 
measurements. 

One  of  the  challenges  of  uantifying  the  turbulence 
intensities,  and  comparing  to  the  measurements,  is  the  fact  that 
it  is  challenging  to  accurately  capt  ire  the  temporal  and  spatial 
variability  represented  by  the  simulation  with  corresponding 
measurements.  At  least  two  factor ;  are  in  play  here.  The  first 
is  the  initiation  and  maintenance  of  the  statistically  stable 
environment.  This  can  be  seen  fr  m  the  variations  shown  in 
Figure  4.  One  may  notice  that  after  initial  setup,  the  variance 
of  both  temperature  and  velocity  >egins  to  reduce  to  a  stable 
level,  with  smaller  oscillations  like  iy  resulting  from  larger  cell 
movements  within  the  tank,  as  rell  as  naturally  occurring 
intensity  fluctuations  local  to  the  sensors. 


Figure  4.  Times  series  of  temperature  and  vertical  velocity  variance, 
taken  from  the  center  line  of  the  tank,  using  cm-resolution  run. 

Obviously,  this  poses  a  challenge  in  quantifying  the 
turbulent  dissipation  rates  we  discussed  earlier,  both  of 


temperature  variance  and  kinetic  energy.  The  energy  spectra 
from  tank  measurements,  and  model  output  are  shown  in 
Figure  5.  We  see  that  they  all  exhibit  a  well  resolved  inertial 
subrange,  where  the  signature  energy  spectrum  slope  of  -5/3  is 
marked.  This  is  a  clear  confirmation  that  we  have  achieved  the 
initial  goal  of  setting  up  a  controlled  turbulence  environment. 
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Figure  5.  Kinetic  energy  spectra  from  lab  measurement  for  both 
strong  and  extreme  turbulence  case  (top),  compared  to  LES  mode 
output  (bottom). 


We  also  notice  from  Figure  5  the  differences  between  mm- 
and  cm-resolution  output,  and  potential  issues  associated  with 
Taylor’s  frozen  turbulence  hypothesis.  More  in  depth 
understanding  and  discussion  is  needed,  along  with  more  case 
studies,  which  will  be  discussed  in  a  separate  paper.  We  now 
turn  our  attention  to  EO  transmission  in  this  controlled 
environment. 


III.  Imaging  Trhough  Controlled  Turbulence 


Figure  6.  Image  transmission  through  non-,  strong-  and  extreme- 
optical  turbulence  from  RB  tank  at  NRL  SSC. 

To  examine  the  variation  of  image  degradation,  a  structure 
similarity  index  metric  (SSIM)  is  used,  which  has  been  shown 
to  be  very  sensitive  to  the  image  quality  variations  with  known 
pristine  reference  [14],  including  time-varying  influences  such 
as  turbulence  impacts.  A  time  series  of  SSIM  from  our  tank 
experiment  is  shown  in  Figure  7.  We  notice  the  reduced 
variance,  which  indirectly  confirms  our  conclusion  of  a 
controlled  environment  of  the  previous  section.  It  can  be  seen 
that  such  metric  does  follow  the  trend  of  the  turbulence 
intensity  reasonably  well,  at  least  for  the  two  cases  of  strong 
and  extreme  conditions  we  tested.  This  further  suggests  that  it 
could  be  used  as  a  proxy  for  OTF,  given  known  optical 
conditions  associated  with  Eq.  2. 


Figure  7.  Time  series  of  image  degradation  during  our  lab 
experiment,  under  strong  and  extreme  turbulence  case. 

In  atmospheric  studies,  where  turbulence  often  affects 
astronomical  observations,  an  effective  approach  is  based  on 
statistical  method  using  layered  medium,  with  diffraction 
propagation  between  layers.  The  phase  variation  of  each  layer 
can  be  easily  obtained  from  the  power  spectrum  of  the  phase, 
in  3d  case,  this  is  typically  [7] 


It  has  been  modeled  and  validated  in  previous  studies  [1,5, 
13],  that  optical  turbulence  does  degrade  image  quality,  as  also 
seen  by  our  tank  experiment  results,  shown  in  Figure  6.  Cases 
of  transmission  through  non-turbulence,  strong  and  extreme 
turbulence  are  displayed.  It  can  be  seen,  along  with  Figure  2, 
that  optical  turbulence  degradation  is  frequency  dependent,  as 
predicted  by  Eq.2. 


tf>(/c)  =  0.033C^  l/cl-11/3  (5) 

with  C„2  denotes  turbulence  intensity;  or  in  practice,  such  as  in 
our  test  case  shown,  modified  von  Karman  spectrum  is  more 
appropriate,  which  eliminates  the  issues  at  very  high  and  low 
spatial  frequencies 


<Pv(k)  =  0.033 


(6) 


where  km  and  k0  corresponds  to  inner  and  outer  scales. 


Slope  of  -5/3 


spectrum  model  of  different  spectral  slopes.  Top:  Kolmogorov  type 
with  -5/3  slope;  Middle:  steep  reduction  with  a  slope  of  -2,  resemble 
small  scale  processes;  Bottom:  gradual  decrease  of  energy  at  a  slope 
of  -1.  All  based  on  von  Karman  shown  in  Eq.6. 

Applying  formulas  derived  for  underwater  conditions  [1], 
phase  screens  using  spectral  method  are  generated,  following 
Lane  et  al  [15].  Typical  underwater  conditions  corresponding 
to  natural  conditions  are  used,  with  inner  and  out  scales  set  at 
0.001  and  1.5m  to  reflect  the  tank  conditions.  The  results  are 
shown  in  Figure  8.  As  mentioned  earlier,  it  is  interesting  to 
observe  the  uncertainty  in  the  slope  of  the  energy  spectrum, 
which  can  span  a  significant  range,  even  in  a  well-established, 
fully  developed  turbulence  environment  (Fig.5).  Coupled  with 


the  progress  and  time  taken  to  reach  a  stationary  environment, 
it  is  reasonable  to  assume  that  different  spectral  slopes  exist, 
and  worth  the  investigation  of  impacts  under  such  variable 
conditions.  This  is  the  reason  behind  the  middle  and  lower 
panels  of  Figure  8,  where  the  influence  of  spectral  slope  is 
shown  in  the  form  of  a  typical  phase  screen,  compared  side  by 
side  to  the  traditional  Kolmogorov  one  shown  in  the  top  panel 
of  Figure  8.  Strong  contrast  can  be  seen,  especially  between 
fast  energy  reduction  (middle)  and  slow  diffusion  (bottom). 
Notice  that  the  color  is  added  to  contrast  variation  and  only  the 
pattern  variability  should  be  taken  into  considerations. 


Figure  9.  Cross-sectional  view  of  temperature  variations  at  the  center 
of  the  tank  under  a  fully  developed  turbulence  condition,  taken  from 
cm-resolution  numerical  results  at  t=850s,  875,  and  950s  respectively, 
from  top  to  bottom. 


To  provide  a  comparison  of  spatial  variability  across  the 
tank  cross-section,  the  temperature  profiles  from  mm- 
resolution  LES  output  are  shown  in  Figure  9,  at  time=850,  875 
and  950s.  One  can  clearly  see  the  large  to  small  scale  vortices. 
Notice  that  the  phase  screen  in  Fig.  8  is  representative  of  the 
center  of  the  tank  conditions,  while  the  cross-sectional  plots 
shown  in  Fig.  9  depict  a  much  larger  area.  Visual  observation 
suggest  that  the  top  two  panels  of  Fig.  8  are  closer  to  reality 
than  the  last  panel,  when  Fig.  9  is  used  as  the  benchmark. 

IV.  Summary 

Efforts  have  been  made  to  create  a  controlled  environment, 
under  which  stochastic  processes  like  ocean  turbulent  mixing 
can  be  studied.  This  paper  presents  a  first  attempt  to  study  the 
impact  of  optical  turbulence  on  EO  signal  transmission 
underwater  under  such  conditions.  A  FES  numerical  model  is 
used  to  simulate  and  monitor  the  mixing  process. 
Measurements  using  AD  Vs  were  taken  to  validate  numerical 
output.  Comparisons  presented  indicate  a  successful  match  has 
been  found  between  the  simulation  and  measurements, 
although  further  exercise  is  pending,  to  obtain  repeatable 
outcome  at  various  grid  sizes  and  turbulence  intensities. 
Optical  transmission  through  the  RB  tank  are  also  recorded. 
Degradation  by  turbulence  is  related  to  image  quality  metric 
and  may  be  used  as  a  proxy  for  turbulence  intensity 
estimation.  Phase  screens  based  on  the  von  Karman  spectrum 
are  shown,  which  are  in  line  with  variations  of  the  numerical 
tank.  Uncertainty  estimation  of  the  power  spectrum  slope  has 
been  discussed,  and  more  efforts  are  necessary  in  order  to 
increase  confidence  in  the  estimation,  and  establish  direct 
comparison  to  a  complex  phase  screen  model  involving 
scattering  and  absorption. 
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